Soft dilepton production and hard thermal loops * 



H. Zaraket 
Laboratoire de Physique Theorique LAPTH, 
BP110, F-74941, Annecy le Vieux Cedex, France 
(February 7, 2008) 



The Hard Thermal Loop expansion, is an attractive theory, 
but it reveals difficulties when one uses it as a perturbative 
scheme. To illustrate this we use the HTL expansion to cal- 
culate the two loop corrections for soft virtual photons. The 
resulting corrections are of the same order as the 1-loop cor- 
rection. 
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I. INTRODUCTION 

The production of soft photons (energy ~ gT, where 
g <C 1 is the strong coupling constant) in a hot quark- 
gluon plasma at equilibrium has been studied by many 
groups, using thermal field theory |ffl-p|, or a semi- 
classical approach as in &-@]. In this talk I will study 
the case of soft virtual photon (with invariant mass 
Q 2 ~ g 2 T 2 ) at two-loop level using the Hard Thermal 
Loop (HTL) expansion. The two loop diagrams are found 
to be of the same order of magnitude as the one loop re- 
sult. This unexpected result solves the problem of the 
absence of bremsstrahlung in the HTL expansion, while 
bremsstrahlung seems to contribute dominantly when us- 
ing the semi-classical approach 



II. MOTIVATIONS AND PROBLEMS 

The starting point will be the bare thermal field theory. 
The imaginary part of the photon polarization which is 
proportional to the invariant dilepton production rate is 
found to be: 



ImH£ ~ e 2 Q 2 (l - 2nj,(|)) 

for q, q a - gT . 



e 2 q 3 T 2 



(1) 



The additional power of g in Imn^ arises from the fact 
that the soft photon is produced by annihilation of a soft 
quark-antiquark pair (each of momentum Q/2). Other 
processes such as q — > jq are not allowed since the quark 
is massless. Having a soft momentum scale, the HTL 



expansion jo| — demands the use of effective vertices and 
propagators. The soft virtual photon production was one 
of the first applications of the HTL resummation scheme 
H : the authors of || have calculated the imaginary part 
of the one loop diagram of Fig. [I]. In Fig. |l| we have 
depicted some of the cutsQ that give relevant physical 
amplitudes. Cuts [a] and [b] of the one loop diagram in 
Fig. [I] can be approximated by : 



Imn^ M (Q)| [a , b] 



T 
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where p is the quark momentum in the loop, where we 
have used explicitly q ~ gT. 






FIG. 1. l-loop contribution to the soft photon production and 
some of the different cuts responsible for different physical pro- 
cesses; in [b] one of the two fermion momenta is space like while 
the other is time like. 

A few remarks are worth saying about this result: 
• The HTL is built to handle soft momentum problems; 
in the hard momentum limit, it tells that we recover the 
bare theory which is a priori what one needs in most 
cases. For example if we take a fermion with momentum 
P, the inverse of its effective propagator is given by the 
sum of the inverse bare propagator (p) and the one-loop 
correction as shown in Fig. |2|: to evaluate this correction 
one should apply the HTL approximation. Thus in the 
one-loop correction in Fig. ^ we use L + P w L, which is 
justified if P is soft and the loop momentum L is hard, 
but when P becomes hard this approximation is no more 
correct, then the HTL estimate of the one-loop correc- 



* Based on work done in collaboration with P. Aurenche, 
F. Gelis and R. Kobes. Talk presented at the TFT'98 confer- 
ence, Regensburg, Germany, 10-14 august 1998. 



J The thermal cutting rules |ij,[14j is a method that one can 
use to evaluate the imaginary part of a Green function. 
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tion of the effective propagator is incomplete in the hard 
momentum limit. 
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FIG. 2. The inverse of the fermion effective propagator. 

• The processes shown in Fig. [l] come from one loop 
corrections to the bare vertices and propagators. On the 

T 

other hand the integral J — ~ \n(l/g) in Eq. |[ indicates 
gT P 

that the quark momentum runs from a soft scale (gT) to a 
hard scale (T), i.e. we need the hard limit of the one loop 
correction which is incomplete in the HTL approximation 
as we mentioned before. 

• The importance of processes coming from cut [a] in 
Fig plcan be understood from the result of Cleymans 
et al J6|-||] , where they calculated bremsstrahlung photon 
production (Fig. 0) using a semi classical approach: their 



result yields ImIP M (<3) 
order as Eq. [I] and Eq. || 



2 g 3 T 2 which is of the same 



FIG. 3. Bremsstrahlung photon production. 

• The bremsstrahlung photon production (Fig. ||) is 
mediated by the gluon exchange, while that of cut [a] of 
Fig. [j] is mediated by the exchange of quark, and it is 
reasonable to think that a complete calculation should 
include both processes. The natural question that arises 
is: how can one reconcile bremsstrahlung with the HTL 
scheme? 



III. 2-LOOP DIAGRAMS 

The solution to the problems mentioned in Sect. |l| 
comes when one looks at 2-loop diagrams in the HTL ex- 
pansion. In two loop diagrams we have more powers of g, 
coming from the qqg vertices, then to compensate addi- 
tional powers of g which may come from the phase space 
we have to have a hard phase space. In the catalogue 
of 2-loop diagrams we skip those which contain vertices 
with no bare analogue: these vertices turn out to be neg- 
ligible when one of the vertex lines becomes hard. The 
remaining diagrams are shown in Fig. ^. As discussed in 
F. Gelis's talk calculating two loop diagrams in the 
HTL expansion , one should pay attention to the coun- 
terterms in order to avoid double-counting (see also |l6| ] 
for more details.). 
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FIG. 4. Two-loop contributions involving bremsstrahlung pro- 
cesses. A black dot denotes an effective propagator or vertex. 
Crosses are HTL counterterms. 

We explore the case where all the quark momenta are 
hard. Assuming hard quark momentum we can replace 
the effective vertices by the bare ones, and take the bare 
propagator for the quark. The simplified version of Fig. ^ 
is represented in Fig. |^. In the same figure we have de- 
picted the relevant cuts that enable one to calculate the 
corresponding contribution in the framework of the ther- 
mal cutting rules. 




R 

FIG. 5. Simplified 
bremsstrahlung processes. 




two-loop contributions involving 



In what follows we will take the photon to be static i.e. 
= (<7o, 0). For virtual photon the production rate per 
unit time and per unit volume is related to the imaginary 
part of the retarded polarization tensor of the photon as 



dN 
dtdx 



dq Q dq a 
' 12tt 4 Q 



— n B (q )lmU 
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The gluon propagator appears linearly in Imn^, 
we can write: 



imn/ 



imn/ 



1 2 <0 + lmn/ | L 2 >0 



(3) 



then 



(4) 



where L is the momentum flow through the gluon. 
Such separation is motivated by the fact that regions 
L 2 < and L 2 > contain different physical processes. 
While L 2 < gives bremsstrahlung of Fig. ||, which is a 
possible candidate to reconcile bremsstrahlung with the 
HTL expansion, L 2 > gives annihilation and Compton 
scattering. Bremsstrahlung is a new process, then we 
can guarantee that it will give a positive contribution 
that should be added to the one- loop result. Moreover 
we can apply different approximations in each of these 
regions. We can treat each region separately and at the 
end we add the two contributions. In the next sections 
we will discuss in detail the region L 2 < 0, while we will 
give some indicative results for L 2 > 0. 
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IV. BREMSSTRAHLUNG 

The matrix element M. of diagrams [a] and [b] of Fig. | 
can be written as (see Jl(| for the complete expression) : 



M ~ M X L 2 +M 2 



(5) 



Writing M. in this form exhibits some of the terms that 
are absent at one loop level in the HTL expansion. In- 
deed the gluon appearing in the effective propagator and 
vertices in the one loop diagram are taken to be light- 
like (L 2 = 0) by construction of the HTL resummation 
scheme, then Mi will be a pure 2-loop contribution. 



A. Extraction of logarithmic behavior 



We can obtain analytically the leading logarithmic be- 
havior in ln(l/g) of Imll^ at 2-loop order. By momen- 
tum power counting we can isolate the term which gives 
the logarithm: it happens that this term comes from A4i, 
enforcing the previous expectation that bremsstrahlung 
contribution is a pure two loop contribution. 
The contribution of this leading term can be written 
schematically as: 



In, 11 ■„"(,,„. 0) = A x e 2 5 2 ( i] m\ f dp j j 



1o 



f(m g ,q ) 



(G) 



Comments: 

• A is a constant that depends on the color factors, and 
it has no g dependence. 

• e 2 comes from the qq^ couplings. 

• g 2 comes from the qqg couplings. 

• comes from the uncut quark-propagator denomi- 
nators, for example the momentum R in Fig. ^| can be 
approximated by: R 2 — 2q a p. 

• The soft gluon thermal mass m 2 = g 2 T 2 [N + N F /2}/9, 
comes from the gluon spectral density. 

• The ultraviolet cutoff T is provided by the statistical 
weights. 

T 

• The integral J dp indicates the hard p behavior. 

o 

• f(m g ,q ) is a function of dimension one. This func- 
tion depends on both q a which appears as a kinematical 
infrared cut-off in the integral over dl and m 2 which ap- 
pears in the gluon spectral density and can also play the 
role of an infrared cut-off for the same integral. 

T 

• The integral J y points out that I vary between 

a soft scale (/ ~ gT) and a hard scale (T). 
The cutoff / can be simplified by taking q <C m g ~ 
gT <C T, then the larger cutoff m g will be the relevant 
cutoff. After this simplification we get: 



T tt« u, ^ WC F e 2 g 2 m 2 g T, ( T 2 

Imn /( 9o ,0) « 5- In — 

07r J q \ m g 

The production rate is then given by: 



(7) 



dN 



dtdx 



brcmss 



dq dq 
8tt 6 



NC C 



: (5>?) 



1o 



In 



<7o 

n 

n> 2 



(8) 



where the sum runs over the flavor of the quarks in the 
loop (e/ is the electric charge of the quark of flavor /, in 
units of the electron electric charge) . 
Numerical estimates of the unapproximated expression of 
Imll^ verify that the condition q a <g; m g is only a tech- 
nical condition as we get a reasonable accuracy even for 
<Zo/ TO g ~ 10- On the other hand m g /T must be smaller 
than 0.1 in order to have an agreement between the ap- 
proximate and the complete expression with an accuracy 
better than 5 % (see Fig. |). 



FIG. 6. Comparison of numerical estimates of the complete 
matrix clement with the simple theoretical expression obtained in 
Eq. (^). Both plots show the ratio "Numerical/Theoretical". On 
the left plot, q /T is fixed at 10 -4 and we look at the variations 
with rrig/T (i.e. with g). On the right plot, m g /T is fixed at 10 -2 
and the photon energy varies from ultra-soft energies to hard ones. 



B. Comparison with other approaches 

1. Braaten et al. results 

The production rate of soft static photons has been 
evaluated by Braaten, Pisarski and Yuan in JH| at one- 
loop in the HTL scheme. To compare Eq. |pwith the 
result of H in the domain q Q -C m g ~ gT <C T for 
which our expression is justified, one has to isolate from 
their result the processes which have the same nature 
and spectrum (in 4r) as the bremsstrahlung, i.e. those of 



B-fal- 
;. EVA 



Fig 

Sect. (IV A 
result 



Applying the same approximation we did in 
we can easily extract the ln(l/g) in BPY's 



dN 



dtdx 



1 — loop 



dq dq 
12tt 4 



Na 2 



(I>?) 



In 
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(9) 



where m 2 F = g 2 C F T 2 /8 is the soft quark thermal mass. 
Comments: 
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The production rate in Eq. || and Eq. || have the 
same spectrum in 

m 2 F plays the role of an infrared regulator in Eq. 
while TO g is the regulator in Eq. This reflects 
the different nature of both results, an exchanged 
fermion in the first, and an exchanged gluon in the 
second. 



• Comparing the two results we obtain the ratio: 



dN 



brcmss 



32 N + N F /2 



dN\ 



3tt 2 



(10) 



1 1— loop 

which for 2 light flavors and 3 colors becomes 



dN\ t . 



loop 



32 

72 



3.2 



(11) 



This ratio is rather large, which means that 
bremsstrahlung is definitely an essential contribu- 
tion to the soft static photon production rate by a 
hot plasma. 



2. Cleymans et al. results 

We return to the semi-classical treatment used in 
In their approach they took into account the effect of 
multiple scattering of the quark in the plasma. To com- 
pare the two loop result with their result one need to 
"undo" the effect of rescattering and consider only one 
collision of the quark in the plasma. We apply the same 
simplification as they did i.e.: 

• The scattering of the particles is treated as in vac- 
uum (the dynamics), while the plasma effects are 
introduced for the kinematics only. 

• The energy of the quarks or gluons is much larger 
than the temperature so that Boltzman distribu- 
tions are used for particles entering the interaction 
region and a factor 1 is assigned to those leaving 
it. As a by-product one can neglect the photon 
momentum compared to the momenta of the con- 
stituents in the plasma. 

• To screen the forward singularity they introduced 
a phenomcnological Debye mass m D , which we will 
take to be m g when we compare their result with 
Eq. g. 

Applying the above approximations we get for the pro- 
duction rate of the lepton pair at rest the following ex- 
pression: 



dN 
dtdx 



dqodq 22,, 
- _ a a„dtd 



3tt 6 



\ " 



nvl. 



(12) 



where df = 2 s x3 c = 6 and d = | 2x2/x2 s x3 c +2 s x8 c = 
26+| arc the degeneneracy factors introduced in Q . 
From Eq. [l2] we see that the semi-classical approach give 
the same spectrum in as expected. 
Comparing with Eq. @ we find for two flavors: 



dN\ 



i semi — c 



lass 



dN\ 



brcmss 



15 

72 



(13) 



i.e. the semi-classical over-estimates the rate of produc- 
tion by about 50%. This difference appears to come from 
the various approximations they used n. 



V. PRELIMINARY RESULTS FOR TIME LIKE 
GLUON 

In this section we give partial results for the L 2 > 
case |l!| and the arguments that we will give are based 
on the leading logarithmic behavior only. 
The basic physical processes for L 2 > are Compton 
and annihilation, then naively we expect double counting 
with the one loop diagram, in the region where the imag- 
inary part of the one loop diagram gives similar processes 
(Fig. 0-[b]). To cure this problem of double counting , we 
have to calculate the counterterms depicted on the right 
of Fig. ||. Indeed, when the gluon in the loop becomes 
hard, we have a hard loop that reproduce what is already 
included in the one-loop diagram via the effective vertices 
and propagators . 

The structure of Imn p M allows to decompose it as the 
sum of two pieces, one part reflecting the hard gluon mo- 
mentum behavior (Imlli) where we can apply the hard L 
approximation, and the other reflecting the hard fermion 
behavior (Imll2): 

Imn ~ Imhi +Imn 2 . (14) 



A. Hard L region 

By isolating hard L terms we apply indirectly the HTL 
approximation, then one may expect a fair compensation 
between Imlli and the counterterms. Imlli reveals soft 
p sensitivity, which means that it is no more allowed to 
use the bare propagator for the quark, thus we have to 
use the effective one. At leading logarithm order we can 
take a constant mass for the quark which will be the 
asymptotic quark thermal mass which is twice the soft 
quark thermal mass m F . This mass will play along with 
go the role of an infrared regulator. 



2 The Boltzman approximation is a favorite candidate of such 
an overestimate. 
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Then the contribution of Im 111 to the production rate is 
approximated by: 



dN 



d 4 xdq a d 3 q 



I nil 



2 2 2 



I2ir 4 



Q 2 



In 



T 



g(m F ,q ) 



(15) 



where g is a function of m/ and q Q . To evaluate the imag- 
inary part of the counter term diagrams (CT) we follow 
fl|, where the authors relate the imaginary part of the 
effective vertex to the spectral density function. Then 
the imaginary part of the counter term diagrams may be 
expressed in terms of the spectral density of the quark 
propagator. After a cumbersome calculation, the loga- 
rithmic behavior of the counter terms can be estimated 
by: 



dN 



d 4 xdq Q d 3 q 



CT 



N T,f e / c?m\ 



12tt 4 



In 



2Tq 



h(m F ,q ) 



(16) 



where h is a function of mp and q a . 
The only difference between Eq. |l5| and Eq. [l6] is that 
they have different infrared regulator, adding both Eqs. 
we get a logarithm of argument of 0(1), which remains 
finite, even in the limit of vanishing thermal masses [], 
then we can neglect the total contributions coming from 
ImlTi and the counterterms compared to ln(l/g) that we 
will get from Ini2 . 



B. Hard p region 

In the hard p limit we are lead to the same type of 
analysis as for L 2 < where there is no need for the 
effective propagator, thus we use the bare propagator for 
the quark. The different behavior of the spectral density 
for L 2 < and L 2 > leads to different behavior of the 
terms in the matrix element, for example the term that 
gives the leading logarithm in Im II2 comes from A^2 in 
Eq. ||. Different behavior will lead to different spectrum 
in each case. 

Taking a constant thermal mass for the gluon, we can 
obtain analytically the leading logarithmic behavior of 
Im 2 : 



Imn 2 (<7o,0) « — — — q T\n 



The production rate is then given by: 




(17) 



A limit to be used later on. 



dN 



d 4 xdq d 3 q 



hard p 2NJ2e 2 f a 2 r 



3tt 6 



In 



T 



( 2k 1 qo_\ 
(18) 

which has a spectrum compared to the of 
bremsstrahlung. 



C. Comparison with Braaten et al. result 

To compare Eq. ^ with the one loop result obtained in 
P|, we consider only the kinematical region which gives 
rise to processes shown in Fig. 0-[b], i.e. the region where 
one of the quarks in the loop is space like while the other 
is time like. Applying the same approximation done to 
obtain Eq. [l7] we get for the diagram [b] of Fig. [I]: 

Hoop 



dN 



d 4 xdq Q d 3 q 



~YhT 4 qj 



In 



( 2Tq 



\Max(q 2 , m 2 



(19) 

which has the same spectrum in as the production rate 

at two loops in Eq. Adding Eqs. [l8| and [l9| we get for 
the annhilation and Compton scattering up to two loops: 



dN 



NJ2e 2 f a 



f 



d 4 xdq d?q 3?r 6 



2Tq 





4 \M&x(q 2 ,m 2 F 



+2 In 



T 



1 go N 

\q ~ m g j 



(20) 



The ratio of the two-loop contribution to the one-loop 
result is given by: 



dN\ 



2 — loops 



dN\,_ 



0.81 . 



(21) 



loop 



This signifies that the two-loop contribution is essential. 
Being calculated with p taken hard from the beginning, 
Im II2 will complete the contribution coming from the 1- 
loop result which assumes p to be soft compared to the 
momentum of the hard loop in the resumed propagators 
and vertices. This cures (at least at leading logarithmic 
behavior) the problem mentioned in Sect. |l[ that the 
hard thermal loop does not give the complete estimate 
of the loop corrections to the effective quantities in the 
limit of hard momentum. To sum-up, the one-loop dia- 
gram takes care of the soft p region, while the two-loop 
diagrams with the counter terms complements the one 
loop estimate by taking care of the hard p region. 

To verify the above result we take the thermal masses 
to zero in the logarithms in Eq. |20|: 



dN 



NJ2ej a 2 g 2 C F( n 2 



d 4 xdq d 3 q - i^^- + 2 ^ T M- ^ 

This coincides with the production rate calculated in Q 
at two-loop level using the bare thermal field theory. 
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VI. FURTHER PROBLEMS 

Our result reveals also a logarithmic sensitivity, i.e. 
the gluon which is supposed to be softQ is extrapolated 

T 

by the logarithmic behavior (J y) to the hard region, 

which bring again the problem of the hard limit of the 
HTL approximation at next to leading order. A complete 
result may need to include higher order contributions as 
the 3-loop diagrams of Fig. fj]. Probably one has also to 
consider the other topologies of 3-loop diagrams to have 
a gauge independent set of diagrams, without forgetting 
the associated countcrtcrms. 



rate in the plasma is increased considerably which may 
have remarkable experimental effects. 




FIG. 7. A set of 3-loop diagrams that may give a dominant con- 
tribution, diagrams where the gluon loop is replaced by a fermion 
loop is also of the same nature. 

The importance multiple scattering of the quark in 
the plasma (Landau-Pomeranchuk effect) in |]-^], points 
again the importance of higher order diagrams. 



VII. CONCLUSION 

We have discussed the production of virtual soft dilep- 
ton production in a hot quark-gluon plasma in thermal 
equilibrium. This special case may give some indicative 
properties of the hard thermal loop expansion as a per- 
turbative scheme. We have shown that to get the total 
contributions at a given order in the coupling constant g, 
one has to add contributions that comes from different 
loop orders, i.e. increasing the number of loops does not 
necessarily mean that the result is suppressed by addi- 
tional powers of g. The case of soft virtual photon we 
considered illustrates that, even in the absence of any di- 
vergences, loop-order mixing occurs. Moreover the two- 
loop diagrams give rise to new physical processes which 
do not appear at one loop. Using the method of coun- 
terterms we give an explicit way to calculate higher loop 
diagrams using the HTL expansion. There exists other 
methods in the literature, such as using a cutoff p0|-p3"| in 
order to obtain the correct behavior of the effective prop- 
agator in the hard momentum limit and to avoid possible 
double counting. Finally adding the contributions coming 
from L 2 > and L 2 < the virtual photon production 
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